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Abstract. Heuristic arguments and numerical simulations support the Belinskii et 
' al (BKL) claim that the approach to the singularity in generic gravitational collapse is 

■ characterized by local Mixmaster dynamics (LMD) . Here, one way to identify LMD in 



o 



collapsing spatially inhomogeneous cosmologies is explored. By writing the metric of 
' one spacetime in the standard variables of another, signatures for LMD may be found. 

Such signatures for the dynamics of spatially homogeneous Mixmaster models in the 
variables of t r (l)-symmetric cosmologies are reviewed. Similar constructions for U(l)- 
■ symmetric spacetimes in terms of the dynamics of generic T 2 -symmetric spacetime are 

0^ , presented. 

o 

m 

O ■ 1. Introduction 

O 

^ Vincent Moncrief and I became interested in spatially inhomogeneous cosmological 

spacetimes as graduate students at the University of Maryland more than 30 years ago. 
Kuchar [1. had extended Misner's minisuperspace approach to quantum cosmology j2] 
to cylindrically symmetric spacetimes containing gravitational waves (Einstein-Rosen 
waves |3J. At the same time, Gowdy had recognized that essentially interchanging 
the radial and time coordinates in the Einstein-Rosen solution yielded cosmological 
spacetimes containing gravitational waves [I]. Moncrief and I had many discussions 
about the Gowdy models — especially those with T 3 spatial topology. The classical and 
quantum mechanical properties of the polarized subclass of Gowdy T 3 models eventually 
became my thesis Even then, we recognized that Gowdy spacetimes are ideal 
"theoretical laboratories" to test formalisms to develop potentially provable conjectures. 
Later work by a number of authors, e.g. |HllZj,has demonstrated the usefulness of these 
models. 

Both Moncrief and I returned to the Gowdy models from time to time, e.g. as 
an example of extendable spacetimes [Hj or quantum field theory in curved spacetime 
Our paths and the Gowdy models crossed again in 1993 at (K)ITP when we used 
generic Gowdy T 3 spacetimes in the collapsing direction to investigate the approach 
to the singularity with numerical simulation ^H|- This project grew to involve others 
(for a summary see [TT] : for a review see [T2]) and explored a variety of collapsing 
spatially inhomogeneous cosmologies. In this Chapter, I will discuss one aspect of my 
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collaboration with Moncrief — the signature of local Mixmaster dynamics (LMD) (see 
for example fT5|). 

The singularity theorems of Penrose, Hawking, and others state that reasonable 
matter evolving from regular, generic initial data will develop singular or otherwise 
pathological behavior if the gravitational field becomes sufficiently strong. The theorems 
do not, however, describe the nature of these inevitable singularities. Penrose's 
cosmic censorship conjectures state that the singularities must be hidden from external 
observers by an event horizon (weak) and/or not be detectable by a timelike observer 
until he/she falls into it (strong). Specific spacetimes are known with many different 
types of singular behavior. Some yield violations of the cosmic censorship conjectures 
but are for various reasons viewed as non-generic. Collapsing spatially inhomogeneous 
cosmological spacetimes provide an arena for exploring the nature of generic singularities 
and for testing strong cosmic censorship. An excellent review of singularities and cosmic 
censorship along with references to the original papers may be found in I15j. 

Long ago, Belinskii, Khalatnikov, and Lifshitz (BKL) [TT)| ITTj had argued that the 
generic singularity was spacelike, local, and oscillatory. This means that, eventually, in 
a generic gravitational collapse, each spatial point will evolve as a separate universe. 
The time evolution at that point will then describe the approach to the singularity of 
the most general spatially homogeneous cosmology — Bianchi Type IX (or VIII) — the 
Mixmaster model ^H]- Mixmaster behavior has been a focus of study for almost 40 years 
(for a review see At any instant, the state of a collapsing Mixmaster universe may 

be described by a single (BKL) parameter u which encodes the anisotropic collapse rates 
of the spacetime. Each epoch of the collapse is a Kasner spacetime characterized by a 
fixed (in time) value of u. As BKL argued, this fixed value cannot be maintained as the 
influence of the spatial scalar curvature in a generic spacetime begins to dominate. The 
change in u is calculated using conservation of momentum in a bounce off the dominant 
term in the curvature potential PBj- It is found that, for 1 < u n < oo, the value of u at 
the n th epoch is given by 

f u n - 1 u n > 2 

Un+1 = I ~4 1< u n < 2 • (1) 

k u n — 1 — n — 

A combination of numerical simulation jTH] and mathematical analysis [201 has 
demonstrated that the asymptotic dynamics approaches arbitrarily close to that 
described by the w-map (P). Note that the sequence {u n } is sensitive to initial conditions 
(due to the subtraction in the denominator). All the epochs with u n — > u n — 1 form an 
era. When u n — > (u n — l)" 1 , a new era begins. The era-to-era evolution of u is described 
by the Gauss map u^ + i = (u^- — [mat]) -1 , where [ ] denotes integer part. The simple 
relationship of the li-map offers, as was known to BKL, an invariant way to characterize 
Mixmaster dynamics at any spatial point (local Mixmaster dynamics (LMD)). BKL 
referred to LMD as oscillatory behavior. 

The other primary type of approach to the singularity in spatially homogeneous 
cosmologies is exemplified by the existence of a final Kasner (or Bianchi Type I) epoch 
(or u- value). This behavior is called asymptotic velocity term dominance (AVTD) since, 
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once in the final Kasner epoch, the spatial scalar curvature never again plays a role in the 
dynamics. The BKL conjecture is that (1) eventually all singularities will be spacelike, 
(2) each spatial point evolves as a separate spatially homogeneous universe, and (3) 
while the most general behavior is LMD, some systems will be AVTD. 

In the past few years, progress has been made in providing a rigorous underpinning 
to the identification of spacetimes with AVTD approaches to the singularity [6,. For 
example, Andersson and Rendall were able to prove the existence of an open set of AVTD 
solutions for cosmological spacetimes on T 3 x R with no symmetries and a scalar field (to 
suppress LMDJ) [22] • Analogous constructions for spatially inhomogeneous cosmological 
spacetimes expected to exhibit LMD do not yet exist. (See, however, Uggla et al [23J). 
So far, the only detailed knowledge of the approach to the singularity in such models 
come from numerical simulations [H]. Numerical simulations of T 3 x R spacetimes 
without symmetry are in progress jS]. These use invariant quantities equivalent to 
a local construction of the w-map to search for LMD. An invariant calculation of a 
quantity equivalent to u will yield a prediction for the next value of u at that spatial 
point. While this approach has been tested in spacetimes of higher symmetry and 
preliminary results show its value in the most general case [34J, I shall focus here on 
another way to characterize LMD that has proven useful in these models |13j . 

Numerical simulations of Gowdy [25:, T 2 -symmetric [2E], an d t/(l)-symmetric 
spacetimes have used (in contrast to the no-symmetry case) variables adapted to 
the system in question. A Hamiltonian whose variation yields the relevant Einstein 
equations may be developed in terms of these variables and their conjugate momenta. 
The variables are chosen so that the Hamiltonian (density), H, has the form of a kinetic 
"energy" plus a potential "energy." Some of the terms in H depend exponentially on the 
configuration variables. These, as well as non-exponential terms, may contain spatial 
derivatives. Neglecting both the exponential and spatial derivative terms typically 
yields a "free particle" Hamiltonian. The solution to the corresponding equations of 
motion would describe the model's velocity term dominated (VTD) limit (if it were to 
exist). A heuristic argument that the model is AVTD would be that substitution of 
the VTD solution into the neglected terms causes them to become exponentially small 
as the singularity is approached. If the VTD solution is not consistent in this way, 
one typically finds two or more of the exponential potentials alternately growing and 
decaying. This yields the (presumably) infinite sequence of bounces that characterizes 
LMD. This heuristic approach, called the Method of Consistent Potentials (MCP), was 
first introduced by Grubisic and Moncrief [2Zj and later generalized 

Following the MCP approach, it was found that apparently different classes of 
potentials were important in the various non-AVTD models. Most strikingly, the 
Mixmaster dynamics involving the standard three minisuperspace potential terms in 
the spatially homogeneous models somehow become bounces off two potentials in £7(1)- 
symmetric collapse and three again in generic T 2 -symmetric collapse. To resolve this 

| It is well-known [23 that Mixmaster dynamics in spatially homogeneous cosmologies will be 
suppressed if a minimally coupled spatially homogeneous scalar field is added. 
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"paradox," one can rewrite the metric of a Bianchi Type IX cosmology in the same 
variables as those of a spatially inhomogeneous model and thus identify the signature 
for Mixmaster dyanamics in these variables [24J. We note that, while we shall focus 
primarily on T 3 spatial topology, the Bianchi Type IX as [/(l)-connection is, of course, 
made for S 3 spatial topology Heuristic evidence is that close to the singularity, as the 
influence of nearby spatial points on the dynamics at any given one becomes negligible, 
topological differences also become unimportant. This may be demonstrated explicitly 
through comparison of the approach to the singularity of Gowdy models on T 3 x R |2S] 
with those on S 2 x S 1 x R gg. 

In this Chapter, I shall summarize the (published) predictions of LMD in U(l)- 
symmetric models and report on some tests of these predictions. I shall then discuss 
the closely related issue of expressing a T 2 -symmetric spacetime as a [/(l)-symmetric 
one. This allows similar predictions to be made from observations in T 2 -symmetric 
spacetimes [28J. We shall see that a number of interesting open questions arise in this 
case. 



2. Mixmaster as a £/(l)-symmetric spacetime 

The material in this section summarizes [T^j. On S 3 x R, one can compare a spatially 
homogeneous (diagonal) Bianchi IX metric as written in a coordinate frame: 

ds 2 IX = - e 2(a+c+7) dr 2 + e 2a (cos dO + sin 9 sin (p dip) 2 

+ e 2C (- sin (p d6 + sin 9 cos (p dip) 2 + e 2l (d<p + cos 9 dip) 2 (2) 

to a spatially inhomogeneous [/(l)-symmetric metric. U(l) symmetric cosmologies on 
S 3 x R are described by the metric [2T?] : 

ds 2 u(1) = e~ 2lp {-N 2 dr 2 + g ab {dx a + iV a dr)(dx b + N b dr)} 

+ e 2lp (dip + cos 9 dcp + f3 a dx a + (3 dr) 2 . (3) 

In (j2J), the logarithmic scale factors (LSFs), a, (, and 7, are functions of r, while (p, 
9, and ip are angles on S 3 . In (jSJ), the symmetry direction is ip, and the other spatial 
directions are {x a } = {9,<p}. The metric variables are assumed to be functions of 9, <p, 
and t. The norm of the Killing field (in the ip direction) is e v , (3 a are the "twists", e 2A 
is the determinant of the 2-metric g a b = e A e a b and e a b is parametrized by x and z via 

e 2z + e~ 2z (l + x) 2 e 2z + e~ 2z {x 2 -l)A 

e 2z + e~ 2z (x 2 -1) e 2z + e- 2z (l-x) 2 \ ' ^ 

Note that the metric differs from that for T 3 spatial topology given in [30 . It 
is convenient to make a canonical transformation from the twists and their conjugate 
momenta e a to the twist potential u and its conjugate momentum r. It is also convenient 
to define the spacetime slicing by zero shift and lapse iVsin# = ^J~g = e A where g is 
the determinant of the 2-metric g a b- In the variables of (J2J , the singularity occurs at 
t = 00. In any Kasner-like, epoch, as r — >• 00, two of the LSFs will be decreasing and 
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one increasing (anisotropic collapse). In terms of the Misner variables [TS] Q, (3±, given 
by 

a = Q-2(3 + , C = Q + (3 + + V3(3_, 7 = ft + /3+ - v 7 ^- , (5) 
the evolution of the spacetime is obtained from the Hamiltonian 

2H IX = o = -pi+pi + P i + v IX (p ± , Q) (6) 

where 

Vix = e 4« + e 4C + e 4 7 _ 2e 2 ^) - 2e 2(c+7) - 2e 2(7+a) (7) 

is proportional to the spatial scalar curvature and P q, p± are cannonically conjugate to 
Q, (3±. Under almost all circumstances, either Vix is exponentially small or one of the 
first three terms on the right hand side of (JJJ) dominates the dynamics. This dominant 
term is associated with the increasing LSF in that Kasner epoch. The increasing LSF 
reaches its maximum value as the corresponding momentum vanishes (and then changes 
sign) — i.e., a bounce off the curvature potential occurs to end the epoch. During this 
bounce, the more slowly decreasing LSF also has vanishing derivative at the bounce and 
then begins to increase. Meanwhile, the most negative LSF changes its time derivative 
at each bounce — continuing to decrease but more slowly. The era ends when this LSF 
starts to increase. 

Rather than the three dominant terms in Vix, the L/(l)-symmetric Hamiltonian 
whose variation yields Einstein's equations is 



Hu(i) = J s3 Hu(i) 



Ne 



-A 



2 , ^- Az 2 , 2 , Aid 2 ^ 2 

%Pz + 2 e Px + gP + 2 e r ~ 2 Pa 



Is 3 sin 9 

+ {(e A e ab ) , ab - (eV) , a A, b +e A [(e~ 2z ) , uX , v -( e ~ 2 *) 
+2eVV,a +^e A e-^e ab u, a cu, b }] (8) 

where H.jj(i) = is the Hamiltonian constraint, the overall trigonometric factor comes 
from N/yfg, and p v , r, p\, p z , and p x are cannonically conjugate to (p, u, A, z, and x. 
The configuration variable u is the "twist" potential obtained from the twists f3 a in the 
metric (fTo^) (see jSOl)- Note the three "potentials" 

V^-rV*, V 2 = ±e A e-^e ab uj, a u, b , V z = \p^. (9) 

The relationship between the Bianchi IX variables a, (, and 7 and the U(l) variables 
(p, 00, A, z, and x may be found by comparison of the corresponding metric coefficients 
in (0) and (jSJ). In this summary, we shall consider only if and z. Figure ^ shows the 
U(l) variables ip and z superposed on the LSFs obtained from a numerical simulation 
of a vacuum, diagonal, Bianchi IX model. Clearly, (p, given by 

e 2 * = e 2a sin 2 9 sin 2 + e 2C sin 2 9 cos 2 + e 27 cos 2 9, (10) 




Figure 1. Evolution of the Mixmaster LSFs a, £, and 7 toward the singularity vs 
t (increasing in collapse). The inset shows the region indicated by the rectangle in 
the upper left hand corner. The U(l) variables ip and z are also shown. Note that ip 
follows the largest LSF while z increases at the end of an era. The arrows show the 
equality of the two subdominant scale factors that marks the era's end. 



tracks the largest LSF. Note that V3 is exponentially small unless z is of order unity 
(rather than large and negative). Since z is given by 

2z 4e 2 («+C+7) sin q ^ e 2 7 cos 2 g + sin 2 g^c cog 2 ^ + e 2« sin 2 ^ 

6 = e 2 ( a +0 + e 2 ( a +T) + e 2 ^) - e 2 ( a +0 cos(2#) + (e 2 ° - e 2 C)e 2 ^ cos[2(# - 0)] ' ^ 

it is clear that, after factoring out the largest (i.e., expanding) LSF, e 2z is of order unity 
when, at the end of an era, the now increasing smallest LSF becomes equal to the middle 
LSF (as seen in Figure 1). This behavior of (p and z, when it appears at a given spatial 
point in a simulation of L^l) -symmetric collapse, may be thus considered a signature 
of LMD. 

Since such U(l) simulations to date are limited in initial data ansatz, spatial 
resolution, and evolution time r, it is not surprising that the signature for the end 
of an era (increasing z) may not appear at all in any given simulation. Figure El shows 
the generically seen oscillatory behavior of (p (as discussed in detail elsewhere |30j). 
The number of (^-bounces shown is typical providing further evidence that the typical 
evolution seen in the simulations represents the midst of an era. FigureElshows a possible 
era-ending ^-bounce from a simulation with different initial data.§ This appearance of 
both ip and z bounces in the U(l) collapse simulations thus provides heuristic evidence 
for LMD in these models. (These U(l) simulations have imposed T 3 , rather than S 3 , 
spatial topology. As was previously mentioned, one would not expect the signature for 
LMD to depend on spatial topology when sufficiently close to the singularity.) 

§ The presence of two z-bounces in figure [3] is actually consistent with the expectations of the it- map 
(|T|). Any era with fractional part of u between .5 and 1 will yield a subsequent era with integer part 1 
and thus a single bounce. This makes single bounce eras quite likely. 
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Figure 2. Typical behavior of y(r) at a typical spatial point. The evolution 
may be characterized as bounces off the U(\) potentials V\ — \r 2 e 4ip and V 2 — 
\ e A e~ 4v e ah uj, a u>b which are also shown. The typical behavior of z (see figure 6 in 
[HOj) is monotonic decrease with decreasing slope where the slope changes when ip 
changes from decreasing to increasing. 
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Figure 3. Evidence for an era end in a U (1) collapse simulation. Note that z increases 
in association with bounces of ip vs r at a given spatial point. 

3. The vacuum Gowdy model on T 3 x R as a £/(l)-symmetric spacetime 

The analysis in this section is based in part on some notes given to me years ago by 
Moncrief. As is apparent from chapters elsewhere in this volume, Gowdy spacetimes 
have proven to provide useful arenas for both numerical and mathematical studies. 
They may be interpreted as cosmologies with spatial topology T 3 , S 3 , or S 2 x S 1 , two 
spatial Killing vectors (KVs), with a natural areal time measuring the change of area 
in the symmetry directions. The gravitational degrees of freedom P and Q depend on 
a spatial variable 9 and time r and represent respectively the + and x polarizations of 
gravitational waves. With 2 KVs, the Gowdy spacetimes on T 3 x R are special cases of 
[/(l)-symmetric spacetimes|| (with 1 spatial KV) on T 3 x R. The relationship between 
the two was developed by Moncrief and used by us as test cases for £7(1) collapse 

| We shall often refer to Z7(l)-symmetric spacetimes as U(l) spacetimes or models. 
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simulations [3T| . 

Einstein's equations for the wave amplitudes P and Q may be found from the 
Hamiltonian jTHI 

H G = \jd6 [nl + e-^vrj + (P 2 +e 2P Q, 2 )] . (12) 
The potentials, 

V, = \e- 2P ^ Q and V 2 = \e~^Q j , (13) 

1 /2 

act to drive the value of v — (ftp + e _2P 7rg) (where P — * vr in the asymptotically 



velocity term dominated (AVTD) limit as r — > oo ^U]) into the range (0, 1) consistent 
with the AVTD behavior of these models. (See [23] for details.) 
Compare the metric in for the Gowdy model: 

ds 2 G = -e x ^/ 2 dT 2 + e x / 2+T / 2 d6 2 + e p -\da + Q d5) 2 + e- p - T dS 2 (14) 

to the metric [32] for the U(l) model on T 3 x R: 

ds 2 u{1) = -N 2 e~ 2v dT 2 + e~ 2 ^e A e ab dx a dx b + e 2 ^{dx 3 + /3 a dx a ) 2 , (15) 

with the gauge condition N = e A . Now identify the corresponding metric components. 

The Killing direction x 3 is identified with the a Killing direction (identification with 
5 cannot be made consistent) so that the norm of the Killing field is 

2 V = P - T . (16) 
The (off-diagonal) U(l) "twists" f3 a are seen to correspond to 

(3 S = Q, Pe = 0, (17) 
while a comparison of g TT yields 2A — 2(p = A/2 — 3r/2 to give 

a-± + *-£. (is) 

4 2 4 v ' 

We require that the conformal metric e a b have unit determinant [32] ■ Since egs = (there 
are no such cross-terms in the metric), we must have ess — e$e ■ The identifications 

e' 2 ^ A e ee = e x ^ 2 and e~^ +A e ss = e~ p ~ T (19) 

yield 

e ee = e A+2 \ ess = e~ A ~ 2r (20) 

so that eee z&8 = 1 as required. This means that the Gowdy slicing is consistent with 
the chosen U(l) gauge condition. 

To construct x and z, we note that the U(l) conformal 2-metric is given by (j4*|). 
Thus, solving for xe~ 2z and e~ 2z in terms of the known components of e a b, we find 

z = -^lncosh(A + 2r), x = tanh(A + 2r). (21) 

We now wish to find the relationship between the U(l) variable u (conjugate 
momentum r) and the Gowdy variable Q (conjugate momentum ttq). Recall that 
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(3s = Q, has conjugate momentum e s = ttq. The U(l) constraint e a , a = is implemented 
by defining the twist potential u through 

e a = e ab u, b (22) 



where 



Therefore 



so that 



e " = (-1 I • < 23 > 



e 5 = cu, e = ttq (24) 



UJ 



d9'TT Q . (25) 



To find the conjugate momentum r, we recall that the kinetic term in the U(l) 
Hamiltonian containing r is^f 

-r 2 e iip and us , T = r e iip . (26) 

Thus 

r = f d& 7c Q , T . (27) 
But the Einstein equations for the Gowdy model yield 

7r Q , T =e- 2T (e 2P Q,e),e (28) 

so that 

r = Q,e- (29) 
The relevant U(l) kinetic term is then 

\r 2 e^= l -{Q, e fe* p -* (30) 

which is the Gowdy potential V<i in (fT3*j) . The £7(1) curvature potential term is 

le A e- 4 ^V,*) 2 - (31) 
But e" = e~ A ~ 2r so that (jSH) becomes 

\-l-- 2P (32) 

which is precisely the Gowdy potential V\ in f)13|) . Note that the canonical 
transformation (f3 a ,e a ) — > (r, u;) has interchanged the roles of V± and V2. 
The remaining momenta are 

p z = Az, T = -2 tanh(A + 2r) (A, T , +2), (33) 

Px = e- 4z x, T =(A, T +2), (34) 
% To obtain the Hamiltonian for T 3 x R C/(l) models, replace sin 8 by 1 in (JHJ ■ 
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and 



so that 



Pa = -A,r • (35) 
The terms in the U(l) Hamiltonian containing these momenta are 

^ = I t anh 2 (A + 2r)(A, r +2) 2 (36) 

^S cosh^W ^ < 37 > 

\p 2 z + \ple 4z = \{^ T +2f. (38) 

Note that if A — > ±oo linearly in r as r — > oo (as expected asymptotically), 
p z — > — 2(A, r +2), a constant, while e 4z — > 0. These behaviors are consistent with the 
MCP analysis of Gowdy spacetimes. 

4. Generic T 2 symmetric spacetimes as [/(l)-symmetic models 

The Gowdy spacetimes are not the most general T 3 x R spacetimes with two spatial 
KVs. As recognized by Gowdy 0], the most general T 2 -symmetric spacetimes have 
additional off-diagonal "twist" metric components. In the vacuum case, the information 
in the twist terms may be reformulated as twist constants (in time and space). Without 
loss of generality, one twist constant may be set equal to zero. The T 2 -symmetric 
spacetimes have been studied analytically jHSj and numerically [2H]. Generic models 
appear heuristically to exhibit LMD with bounces off generalizations of the Gowdy 
potentials V\ and V 2 and off a new twist potential V3. Einstein's equations may be 
found from the variation of the Hamiltonian density 



H T 2 = - — \ ti p + e 7T Q + e [P, 6 +e Q, 



4tt a 

+ ™ 2 7r Ae ( A+2P + 3 ^ 2 . (39) 

In addition, there is a nontrivial momentum constraint: 

Tip P, e +tc q Q, e +tt a X, e = 0. (40) 

Details may be found in [28J. T 2 -symmetric spacetimes are a special case of U(l)- 
symmetric spacetimes. 

The metric for the T 2 -symmetric spacetimes is given by |2H] 

ds 2 T2 = - e^ 3T ^ 2 dr 2 + e^+^W + ae p - T [dx + Qdy 
+ (6?i + QG 2 )d6 - (Mi + QM 2 )e- T dr] 2 

+ ae- p ~ T {dy + G 2 d6 - M 2 e- T dr) 2 (41) 

where the Gowdy metric (|14|) is recovered if the shifts, M a , and twists, G a , vanish and 
fi and a are set = 1. The G a 's and M a 's are related to the twist constant k through 

e~ T M ue +G lyr = KQe^ +2P+ ^' 2 e^\ (42) 
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e~ T M 2 , +G 2 , T = - K ^ x+2P+ZT ^ 2 e^\ (43) 

This metric may be simplified (without loss of generality) by following a procedure 
introduced by Weaver [35: to eliminate the shifts. It is also convenient to revert to a 
and 5 as coordinates. To avoid confusion with the coordinate of the same name (and 
because this is done in our computer simulations), set a = 1 in (}4~Tj) . Let 

a = x- J' air' M 1 e~ T \ 5 = y - f air' M 2 e~ T \ (44) 

and define 

9 = Ke^ x+2P+3 ^ 2 e^ 4 . (45) 

If these transformations are substituted into (j41j) . we obtain the same results as if we 
demand that the shifts vanish — i.e. the M a = in ()42|) and ()43j) . The latter condition 
implies 

d = J" dr' (QQ), G 2 = - f dr' (46) 
so that (|4~Tj) becomes 

ds 2 T2 = - e^-^dr 2 + e^ +T ^ 2 d6 2 

+ e p ~ T Ida + Qd5 + J" dr' {QQ) ~ Q f dr' 6 deX 
+ e- p ~ T [d8-(J T dT'Q)d6} 2 . (47) 
From (}4Tj) . we can identify the U(l) twists as 

fa = Q, (48) 
(3 e = f dr' (QQ) — Q J T dr' 6. (49) 

Note that as r — > oo, j3g — > since Q becomes constant in r. For future reference, we 
note that 

Ps,t = Q,t, Pe, T =-Q,T fdr'Q. (50) 

We shall see below that these assignments are consistent with the relation between the 
f/(l) twists and ui given most simply by [3~2] 

(3 ajT =-Ne-^e ab e bc u;, c . (51) 

As in the Gowdy case, we identify (see ()24|) ) u,g = ttq. 

Now consider the 2-metric g a b which we define here to include all factors in (}4Tj) : 

m = e (A+M+r)/2 + ( r e y e -P- T) (52) 



gss = e p T , (53) 
9se = -e- p - T Te. (54) 
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Note that 

9eo988-9l = ^ T) ' 2 ^e^ (55) 

so that (jnnj) defines A. 

Since the norm of the Killing field in the a direction is g aa = e p ~ T , we identify 

2<p = P - r (56) 

as in the Gowdy case. Combining (f53J) and (JHHj) yields 

A P 5r a ._. 

A = - H h -. (57) 

4 2 4 4 v ; 

Note that the identification of A is different from that in the Gowdy case. 

We also note that this implies that the U(l) gauge condition used previously and 

in [SOI must be generalized. A comparison of (fTKj) and (jUj) yields 

jV = e (A + 2P-5.)/4 = e - M /4 e A (5g) 

so that N ^ e A which violates the previously used U(l) gauge condition. In the 
following, (JH%|) shall be used as the gauge condition. 

The generalized T 2 symmetric models satisfy [2E| = ^e^ 4 so that N = e A /(27r A ). 
This choice of lapse will provide the needed 7Ta's in the denominator of the Hamiltonian 
density (jHUj) and strongly suggests that a similar gauge condition would be needed to 
describe any generalization of Gowdy spacetimes (e.g. magnetic Gowdy JSHI ) as U(l) 
models. 

Now identify the components of the conformal metric e a b using g ab = e A_2</3 e a b to 

find 

e ee = e (A+2P + 3r +M )/4 + e -(A + 2P+3r +M )/4 ^ ^ Q V ^ (5g) 
ess = e -(A+2P+3r +At )/4 ; (6Q) 



e 5fj = _ e -(A+2P+3r +M )/4 / dT , Q (61) 

We shall also need e ee = e -( A + 2P + 3r +^/ 4 . Using fllj) and (EZJ) with CT-CT gives 

e -2z = e A+2r + e -A-2r A + T ^ @\ ' (g2) 



and 

e A+2r _ e -A-2r []_ _ ( f T @^2" 

(63) 



X 



e A+ 2r + e -A-2r( 1 + J- Q )2 



Now consider the Q degree of freedom. From (|51)L we require 

/3e,r= -Ne^eeg^u^) and /3 J)T = -Ne'^eggcu^ ■ (64) 

Our previous identifications for the £7(1) twists in (|49|) and (J4*%|) and our identifications 
of the other variables yield the required relations 

Q ;r = e - 2P ^ and -Q, T f T dr'Q = -e~ 2P ^ f T alr'Q. (65) 
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The term containing r in the Hamiltonian becomes (with our new lapse condition) 



N 1 r 2 

e A 2 4tt a V ' 

so that the variation yields 

r = 2ir x e- Av J e dT'ir Q , T . (67) 
But, in these models, 

^ =e de\-^r) (68) 

so that r = Q,g as in the Gowdy case and 

2 p 2(P-r) 

^e 4 ^ = - (Q, e ) 2 . (69) 

4tt a 4tt a 

The other Gowdy-like potential term is 

(Ne A )e~ A e ee e- Av (uj, e ) 2 . (70) 

Following the steps as with the Gowdy model (since e ee is the same as in that case) 
yields for this term 

^e~ 2P . (71) 
4vr A 

Since the two terms in the Hamiltonian containing Q arise here essentially as in the 
Gowdy case, we must look elsewhere for the twist contribution. 
Let us construct 

;Pl + 7;Py z ) (72) 



e A \8 rz ' T 

from the £7(1) Hamiltonian. Note that our gauge condition means that N/e A = e _/i//4 = 
l/(27r A ). Since z and x are known, the variation of F yields expressions for p z and p x 
in terms of z, T and x, T . We find 

2A+4r 



p x = — I e -2A-4r +M /4 



2 J T dr'Q-(J T dr'Q) 2 ^ 6 



V 



-l+e z 

-2e 2A+4r (l + f dr' e) (A, r +2) 
2e^ 4 [2 (1 + f dr' 6) 6 + (-1 + e 2A - 2 / T dr' 6 - (/ T dr' 0) 2 ) A, 



(73) 



(74) 



(-1 + e 2A + 2 J T dr' + (/ T rfr' 6) 2 ) 
Remarkably, while neither term on the right hand side of ()72|) looks simple, we find that 

F = I e - 2A - 4 ^/ 4 6 2 + ie^ 4 (A, T +2) 2 . (75) 

The second term on the right hand side of (ffH|) will reproduce the appropriate 
generalization of the comparable Gowdy term (j36J) . On the other hand, substitution 
of (|43j) . (jHHj) . and the definition of /1 in terms of 7r A gives 

l e -2A + 4r +M /4 Q 2 = 7TAft .2 e (A + 2P+3r)/2 (?6) 

which is precisely the twist potential. This means that twist bounces (i.e. bounces off 
the twist potential) are associated with ^-bounces in these models. 
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Figure 4. Generic T 2 -symmetric collapse simulation in terms of U(l) variables ip 
and z at a representative spatial point. The arrows mark the twist bounces (see |28) 
for details.) Equation l|56|) defines ip in terms of the T 2 -symmetric variable P and 
essentially follows its behavior. Note that z starts to increase at a twist bounce as well 
as at the end of an era. 



5. Discussion 

In Section 2, we reviewed the analogies between Bianchi IX and [/(l)-symmetric models 
to argue that "^-bounces" should occur in U(l) models at any spatial point that 
experiences the end of a BKL era. Numerical evidence was presented in figure El Here 
we look for the analog of twist bounces [2H] in L^l) -symmetric models. Figure |U shows 
if and z (as functions of P, Q, A, and r) from a T 2 simulation with the twist bounce 
marked. Thus the signature for a twist bounce in a U(l) model should be a change in 
the slope of <p accompanied by a change in z from decreasing to increasing. 

To explore U(l) simulations for these features, we must first consider the gauge 
condition. To understand the difference between the our original gauge condition, 
iV = e A 132 and our revised gauge condition N = e A /(2TC\) ([58)1 . we digress to the 
homogeneous cosmology case. Recall that for anisotropic LSFs a, (, 7, we define 
3f2 = a + £ + 7. Einstein's equation for Q is 



Here we distinguish the 3-D lapse M from the 2-D lapse N used in discussion of U(l) 
models. In the homogeneous models, the BKL time r is defined by the gauge condition 
jV = y/^g. This is equivalent to the original gauge condition since, from the U(l) metric 
(I15j) . Af = e^ v N while \f^g = e~ v e A so that the 3-D gauge condition is maintained. 

On the other hand, if we wish to use a geometrical (e.g. volume) time coordinate, 
then (e.g.) we may choose Cl — — 1 which is equivalent to the gauge condition 



Af 



(77) 




Af = e-v/pn. 

If pn is a constant, the two gauge conditions are equivalent. 



(78) 
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Figure 5. A possible twist bounce in a U(l) collapse simulation using the modified 
gauge condition N = e /pa- The U(l) variables (p and z vs r are shown at a single 
spatial point. The inset displays the feature marked by the rectangle and shows a 
change in the slope of (p along with a change in z from decreasing to increasing, the 
signature for LMD in T 2 -symmetric collapse. The vertical lines are used to indicate 
the alignment of features in tp and z. 

The Gowdy areal time choice is then comparable to the volumetric gauge condition 
fl78j) . (Note that in jS] it is shown that writing (polarized) Gowdy models as 
inhomogeneous generalizations of Bianchi I cosmologies in Misner's minisuperspace 
variables requires a rotation such that — p^ + p 2 + — ► —p T P\ so that the Gowdy spacetime 
equivalent of (J77J) will be f = (A/"/ \f^g 7r A .) Since 7r A is constant in the Gowdy spacetime, 
the areal time choice is compatible with the BKL gauge condition. Once we consider 
generalized T 2 -spacetimes, the BKL gauge condition is no longer compatible and we 
must use the areal gauge condition which is the appropriate analog of (|75|l. 

To better search for twist bounces in U(l) models, we replace the gauge condition 
N = e A with N = e A /p A as an analog of (|78|). Figure El shows a possible twist bounce 
in such a simulation although one cannot regard the evidence to be compelling. The 
required signature is present although at a barely detectable scale. 

Nonetheless, we may argue that this approach of defining signatures of LMD in 
simpler models has proven (and will prove) useful. Future work consists of more careful 
analysis of improved U(l) simulations as well as the development of a direct connection 
between non-diagonal Bianchi IX [37] and T 2 -symmetric models. Discussion in |28j 
indicates that the LMD seen in T 2 -symmetric models can be understood as requiring a 
bounce off a rotational potential wall. Such a wall requires generalization of the metric 
in (J2J) and ^Hj to include non-diagonal metric components for the metric expressed in 
spatial 1-forms. 

The need for the new gauge condition and the need for a more general Mixmaster 
model may imply that the U(l) metric used in [HUj should be generalized. The search 
for LMD signatures can provide heuristic evidence for or against the generality of these 
simulations. 
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